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ABSTRACT 



OUDt£Y KNOX UBRARY 
WVAl POSTGRADUATE SCHOOL 
MONTEREY. CA 93940 



A computer data reduction system has been reprogrammed to 
sample individual scan recordings of the line- spread function 
of laser light propagating through the turbulent marine boud- 
ary layer. 

The distribution of spatial frequencies of the Modulation 
Transfer Functions from random single scan samples from 
experimental data sets has been investigated. A preliminary 
analysis indicates that the zero spatial frequency component 
of the MTF is log-normal distributed. 
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I. INTRODUCTION 



Atmospheric turbulence adversely affects the propagation 
of a laser beam or the performance of an imaging system. The 
Modulation Transfer Function (MTF) is a measure of the degra- 
dation of resolution and thus the relative ability of the 
atmosphere to transmit spatial frequency information. Knowl- 
edge of the MTF of the atmosphere is necessary for prediction 
of the beam cross-sectional profile on target and spot wander 
for a high energy laser and its beam projection system. 

Members of the faculty of the Naval Postgraduate School 
(NPS) , Department of Physics and Chemistry, have investigated 
the effects of atmospheric turbulence on the propagation of 
optical radiation in the marine boundary layer. 

The work started in 1973 and has continued through 1978. 
The research has been concentrated on developing a model for 
the MTF depending on wavelength, range and a properly path- 
position weighted value of the turbulence structure constant 
for optical index Cn. The' work for the period 1 July 1974 - 
30 September 1977 is described in Ref. 1. Field experiments 
were performed by transmitting laser light of different wave- 
lengths along various seapaths . The received intensity was 
detected by a slit-scanning telescope system and recorded on 
tape. The Cn value was determined by computer calculations. 

The objective of this report is to describe the develop- 
ment of a computer program which samples individual scan 
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recordings of the line-spread function and determines a pre- 
liminary statistical distribution of the spatial frequency 
content of the Modulation Transfer Function. 
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II. ANALYSIS OF DATA 



A. FUNCTIONAL DESCRIPTION 

1. Data 

Experimental data are the intensities of the line- 
spread functions as recorded on a 4-channel magnetic tape. 

One track is used for verbal identification of date, time 
and pertinent experimental information. Two tracks are used 
for recording scans and the last track is used for recording 
the trigger signal. 

Verbal information has been cross-checked against 
experiment log book and MTF data output. Data have been 
recorded for one and two wavelengths simultaneously. Only 
data for .6328 microns have been analyzed in this report. 

The line-spread function was sampled by the scanning detector 
at a rate of 50Hz and recorded at a tape speed of 37.5 inches 
per second. 

2 . Data Reduction 

The data reduction is done on a Nicolet Instrument 
Company, NIC-80 data processor. This is a data digitizer 
followed by a 12K 20-bit digital computer and a Tektronix CRT 
display unit. The NIC-80 computer is controlled by a Hewlett- 
Packard HP-9825 desk top type programmable calculator. 

The NIC-80 computer reads the line-spread function 
for one scan, performs a Fourier Transform (FT) analysis of 
the signal and transfers the FT spatial frequency components 
0 - 19 to the HP-9825. The MTF of the atmosphere is calculated 
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by dividing the FT of the line-spread function by the FT of 
the optics. The FT of the optics is stored by the NIC-80 and 
transfered to HP- 9825. The experimental analog signal is 
digitized at 1024 points. The DC component is removed by 
manual setting of signal bias based upon observation of scan 
on CRT display. Any remaining DC component is compensated by 
subtracting the average intensity value of 32 points at each 
end of the digitized signal. 

When not performing FT and transfer operations, the 
NIC- 80 remains in signal display mode with successive scans 
being displayed on CRT. 

The computation and transfer time for one scan has 
been observed to be 5 seconds, and will be the time between 
sampling of scans. During data reduction, the tape recorder 
was operated at 3.75 IPS, corresponding to 5 scans per second. 
This implies that about every 25th scan of the recorded data 
was sampled for analysis. The number of scans sampled was 
limited by the length of experimental data recorded, normally 
in the order of one minute, or 3000 scans. This allowed a 
sample size of 100 - 125 for data analysis. An attempt to 
run the recorder at .375 IPS speed was not successful due to 
problems with the trigger signal. 

3 . Statistical Analysis 

The HP-9825 calculates the MTF at 20 different spatial 
frequencies from 0 to 38000 per radian in steps of 2000 per 
radian. These spatial frequency points are labeled 0 to 19 
and are tested for normal and log-normal distribution. Goodness 
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of fit is determined by the Chi Square test. A data summary 
containing mean, standard deviation and Chi Square values 
for normal and log-normal tests is outputted on the printer. 
Mean and + 1 standard deviation for frequency points 0-19 
is plotted, with the zero spatial frequency normalized to a 
value of unity and all other data adjusted correspondingly. 

A histogram showing relative distribution of MTF values for 
one specific spatial frequency is plotted as requested in 
program input. 

Examples of the above are given in Figures 1, 2 and 3 

respectively. The relationship of the MTF distribution to 

the Cn value has been investigated. The Cn values ranged in 

- 8 - 7 

the order of 8x10 to 1x10 . Typical line-spread functions 
for two different Cn values are shown in Figure 4. Note that 
Cn squared values are given. 

The value of Cn used has been taken from existing MTF 
data outputs. If no output was available, an MTF analysis 
was performed. A detailed discussion of MTF and Cn is given 
in Ref. 1. The scan presentation is a function of time and 
detector position. The dwell time setting was obtained by 
analyzing calibration data recorded on tape. Dwell time is 
the sampling time per data point. Data have been obtained 
by recording the detector output with a diffraction grating 
of known dimensions placed in front of the telescope. This 
provides a known reference of angular spread. The dwell time 
is adjusted so that the angular spread of the first order 
diffraction as observed on the display is equal to the ref- 
erence value . 
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B. PROGRAM THEORETICAL CONSIDERATIONS 

The data were tested for normal and In normal distribution. 
1. Normal Curve 

The normal or Gaussian curve is expressed as 
Y = constant exp )^) 

where 

x is the value of MTF for normal test distribution 
or value of In MTF for log-normal distribution. 

x is the unbiased estimate of the Mean of the MTF 
values or Mean of In MTF values . 

s is the unbiased estimate of the standard deviation 
of the distribution. 

Y is the probability density of the distribution. 

Constant = i/ s/2tt , where i is the width of a histogram 
bin to make area under curve equal to i. 

The Mean x is the arithmetic Mean given by: 

N 

x = E x. /N 
i=l 1 

where N is the number of observations (samples) . 

The standard deviation s is the positive root of the 
Variance s squared given by 

2 N - 2 

s Z = Z (x. -x) / (N-l) 

i=l 1 

For computational purposes the following algebraically 

equivalent formula has been used: 

9 N 9 N 9 

s = ( Z (x. ) - ( E x. ) Z /N) / (N-l) 
i=l 1 i=l 1 
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The probability density curve has been computed and 



plotted over a range of Mean + 3 standard deviations. 

2 . Goodness of Fit 

Test of goodness of fit to tested distribution has 
been done by histogram and Chi Square calculation. 

The width of the histogram bins is a program variable 
expressed as a ratio of the standard deviation. 1/3 has been 
used corresponding to a total of 18 bins . The number of 
samples with values within a bin is expressed as a fraction 
of 1, so that the area of the histogram is equal to the width 
of the bin interval. 

The Chi Square goodness of fit value has been calcu- 
lated from 

.. - ? [Nf (n) -F (n) J 2 

x - RFT55 



where 

N = number of observations 
n = number of histogram bins 
F(n) = number of observations within a bin 

The term Nf(n) is the theoretical number of observations 
which should fall within a bin. The probability area has been 
calculated as the sum of a rectangle and a triangle approxi- 
mating the area under the normal curve. The sampling 
distribution of the Chi Square statistics is tabulated as a 
function of degrees of freedom and percentile. The number of 
degrees of freedom is calculated from 
df = k-n-1 
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where 

k = number of bins in histogram 

n = number of parameters estimated 
For a 18 bin histogram with mean and standard deviation 
estimated, the degrees of freedom is 15, with a 95 percentile 
value tabulated as 25 . A calculated Chi Square value less 
than 25 will, at the 5% significance level, imply that the 
data fits the given distribution. 

C. COMPUTER PROGRAMS 

1 . NIC-80 Data Reduction Program 

This program has been developed utilizing subroutines 
from the existing MTF analysis program. Pages 3 and 4 of the 
MTF program have not been changed. Main program is on page 1 
and a data transfer subroutine has been inserted on page 2. 
Program flowchart is given in Figure 5A. Program listing is 
on page 122. 

2 . Data Analysis Program 

Due to limited HP-9825 memory capacity, this program 
has been split into two parts stored on separate program 
files. During operation, programs are loaded and executed 
under program control. Program "A" is on file 8 and contains 
NIC-80 control and all analysis functions except data summary, 
which is on program "B" located on file 9. Programs are 
documented on flowcharts (Figures 6 and 7) . Program listings 
are on pages 117-121. Execution of analysis program requires 
HP tape inserted with track 1 available for data storage. 
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Each data file must be minimum 160 bytes long and the number 
of files must be equal to or larger than the number of scans 
sampled. 

Execution time for a 100 scan analysis run is of the 
order of 90 minutes. Most of this time is required for tape 
movements and operations . 
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Figure 5A. Program Organization 
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Figure 5B. Program Organization 
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Date: 6 figure #, date, month, year 
Time: 6 figure # 

Cn: As argument and power 
Line 18 



For MTF spatial frequencies 0-19 
Line 19-21 



To verify absolute and relative 

values 

Line 22 





NO, or already LN data on file 



Figure 6A. Data Analysis Program 



24 




I 



Label 

Plot 



Relative distribution of MTF 
Line 32-42 



T 



Calculates and plots normal curve 
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Line 43-46 



Plot 

Normal 

Curve 



V 




Figure 6B. Data Analysis Program 
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Figure 6C. Data Analysis Program 
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Line 55-57 



NO 




26 



Line 58 



Line 59 



Line 60-72 



Line 74 



Figure 6D. Data Analysis Program 
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Return point next spatial 
frequency. Line 1 




Figure 7A. Data Summary Program 
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Figure 8B. Data Summary Program 
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III. RESULTS 



A total of 24 data analysis runs have been performed. 
Five of these have been on the same data sets in order to 
investigate consistency of distribution within the same data 
set . 

An arbitrary grading system for goodness of fit over a 
range of spatial frequencies has been established as follows 
All components in range within 95 percentile = 1. 

One component in range outside 95 percentile = 2. 

Two or more components in range outside 95 percentile 

= 3. 

Data summary and MTF mean and standard deviation plots 
are given on pages 39-86, arranged in order of decreasing 
order of Cn. A distribution analysis summary is tabulated 
on pages 31-33. A set of histograms for all spatial frequen 
cies for one data set is given on pages 87~112- 

A. CONSISTENCY OF DISTRIBUTION 

Data set number 12 has been run 5 times to verify consis 
tency of distribution as only every 25th scan was sampled. 

Utilizing the grading system for spatial frequency range 
0, 1-6, 7-12 and 13-19, the following results as given in 
Table II and III, were found: 
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TABLE I 



DISTIRUBTION ANALYSIS SUMMARY 



Set 

No. 


Cn 


Date 


Time 


Spat . 
Freq . 


Norm 


Log 

Norm 


Comments 




1.04 


031176 


2320 


0 


1 


1 


0 freq. log 










1-6 


1 


1 


data best fit 


1 


IQ" 7 






7-12 


1 


1 










13-19 


2 


3 






1.06 


021877 


1225 


0 


1 


1 


0 freq. normal 


9 








1-6 


1 


1 


data best fit 




IQ' 7 






7-12 


2 


3 










13-19 


3 


3 






1.06 


040276 


1420 


0 


2 


2 


Large standard 


3 








1-6 


3 


3 


deviation 




ID' 7 






7-12 


3 


3 










13-19 


3 


3 






1.07 


031176 


2247 


0 


1 


1 


Good fit at 


4 








1-6 


1 


1 


low frequen- 




1 — 1 
O 

1 






7-12 


1 


1 


cies 








13-19 


3 


3 






1.67 


031176 


2205 


0 


1 


4 

1 


0 freq. log 




• 






1-6 


1 


2 


data best fit 


D 


10" 7 






7-12 


1 


2 










13-19 


2 


3 






1.75 


040276 


1203 


0 


1 


1 


0 freq. log 


6 








1-6 


3 


3 


data best fit 




IQ' 7 






7-12 


3 


3 










13-19 


3 


3 






2.79 


092177 


1555 


0 


2 


1 




7 








1-6 


3 


3 






10' 8 






7-12 


3 


1 










13-19 


3 


1 






4.25 


121776 


1140 


0 


2 


1 


log data best 


3 








1-6 


i 3 


1 


fit 


O 


oo 

1 

o 

t — 1 






7-12 


3 


3 










13-19 


3 


3 





Legend: 1. 

2 . 

3. 



All values within Chi Square 95 percentile. 
One value above Chi Square 95 percentile. 
Two or more values above Chi Square 
95 percentile. 
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DISTRIBUTION ANALYSIS SUMMARY 



Set 
No . 


Cn 


Date 


Time 


Spat . 
Freq . 


Norm 


Log 

Norm 


Comments 


5.03 


092377 


1452 1 


0 


1 


2 


0 freq . log 


9 


1 — 1 
o 

00 






1-6 


2 


1 


data, chi sq. 








7-12 


3 


3 


value 26.6 








13-19 


3 


2 






5.62 


021877 


1103 


0 


1 


1 


log data low 


10 








1-6 


1 


1 


chi sq . values 


o 

00 






7-12 


2 


1 








13-19 


3 


2 






5.86 


121776 ! 


1244 


0 


2 


2 


Absolute values 


11 








1-6 


2 


1 


low 


00 

1 

o 
1 — 1 






7-12 


2 


3 










13-19 


3 


3 






5.89 


092377 


1522 


0 


2 


1 


This data set 


12 


M 

O 

1 

00 






1-6 


3 


1 


run 5 times to 






7-12 


3 


1 


check consis- 








13-19 


3 


1 


tency within 
data set 






5.94 


092277 


1616 


0 


1 


2 




13 








1-6 


3 


3 




1 — 1 
o 

00 






7-12 


3 


2 










13-19 


3 


3 






6.08 


092277 


1705 


0 


1 


1 


log data range 


14 








1-6 


1 


1 


0-8 low chi sq. 


1 — 1 
o 

00 






7-12 


1 


2 


values 








13-19 


2 


2 






6.92 


121776 


1012 


0 


1 


1 


log data good 


15 


00 

1 

o 
1 — 1 






1-6 


1 


1 


fit 






7-12 


3 


1 










13-19 


3 


1 






6.99 


021577 


1700 


0 


2 


2 




16 








1-6 


2 


3 




O 

1 

00 






7-12 


3 


3 










13-19 


1 


3 





Legend: 1. 

2 . 

3. 



All values within Chi Square 95 percentile. 
One value above Chi Square 95 percentile. 
Two or more values above Chi Square 
95 percentile. 
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DISTRIBUTION ANALYSIS SUMMARY 



Set 

No. 


Cn 


Date 


Time 


Spat . 
Freq . 


Norm 


Log 

Norm 


Comments 




7.28 


021777 


1900 


0 


1 


1 




17 








1-6 


1 


1 






10 -8 






7-12 


3 


3 










13-19 


3 


2 






7.35 


021877 ' 


1035 


0 


1 


2 


Large standard 


18 








1-6 


3 


1 


deviations . 


00 

1 

o 
* — 1 






7-12 


3 


3 










13-19 


3 


1 






8.23 


021577 


1635 


0 


2 


2 


Large standard 


19 








1-6 


2 


2 


deviations . 




00 

1 

o ! 

» — i 






7-12 


3 


3 










13-19 


3 


3 






8.56 


021577 


1600 


0 


2 


1 




20 








1-6 


3 


1 






00 

1 

o 
» — 1 






7-12 


3 


3 










13-19 


3 


3 





Legend : 



1. All values within Chi Square 95 percentile. 

2. One value above Chi Square 95 percentile. 

3 . Two or more values above Chi Square 
95 percentile. 
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TABLE II 



Grading Consistency of Data from Same Data Set 



Frequency Range 
0 

1-6 
7-12 
13 - 19 



Grade 



Number of Runs with Grading 
Normal Log Normal 



12 3 
14 0 
10 4 
0 0 5 
0 0 5 



12 3 

4 10 

5 0 0 

3 11 

3 2 0 



The number of samples which met the 95 percentile criterion 
for each spatial frequency in the 0-6 frequency range is 
listed in Table III. 



TABLE III 

Consistency of Individual Frequency Components 
from Samples of Same Data Set 



Frequency 

0 

1 

2 

3 

4 

5 

6 



Number of Samples within 95% 
Normal Log Normal 



1 

1 

2 

3 

3 

2 

3 



4 

5 
5 
5 
5 
5 
5 
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The results indicate that the distribution can be expected to 
give same Chi Square values for different sample series from 
same data set. 

B. FACTORS AFFECTING DISTRIBUTIONS 

1 . Scintillation 

If the intensity of a laser beam which has propagated 
through a turbulent medium is observed at a point, the inten- 
sity appears to fluctuate or scintillate in a random way. 

It has been shown that the values of the logarithm 
of the amplitude follow a Gaussian distribution (Ref. 1, 
page 15) . Spatial frequency component zero represents the 
area under the curve (signal) . This area is proportional to 
the received intensity. It is thus expected that the distri- 
bution of the zero spatial frequency component should fit 
the log normal distribution. 

2 . Signal to Noise Ratio 

The amplitude of the spatial frequencies decreases 
with increasing frequency number. For spatial frequency 15 
and above, the amplitude is of the order of 10% of the zero 
spatial frequency component. At higher frequencies noise 
will contribute more relative to the amplitude. Noise is 
considered normal distributed, and it is expected that the 
distribution at higher frequencies is normal . 

C. RESULTS AND DISCUSSION 

1 . Log-Normal Distribution of 0 Spatial Frequency 

In 13 of 20 runs the spatial frequency zero fits the 
Ln normal Chi Square 95 percentile criterion. For runs which 
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do not meet the distribution criterion, the ratio of mean to 
standard deviation is in the order of 2. This may indicate 
saturation, and thus not follow log normal distribution (Ref. 
1, page 16) . 

The results do indicate that the zero spatial freq- 
uency is log normal distributed. 

2 . Normal Distribution of Higher Frequencies 

For data analysis the modulus, e.g., absolute value 
of the data has been used, and sign of the data has not been 
included. This will introduce a skewness of distribution. 

It is believed that this is the reason for lack of normal 
distribution for data at higher frequencies . 

3 . Distribution versus Cn Value 

Data have been arbitrarily grouped into ranges with 
respect to Cn values as follows: 

Cn High: Cn > 6.10" 8 

Cn Medium.- Cn between 2.10 8 and 6.10 8 

Cn Low: Cn < 2.10 8 

The ratio of distributions graded 1 for the various 
Cn ranges is listed in Table IV. 

The results do not indicate that the distribution 
does change with Cn values. For all analysis results, it 
should be noted that the sample size is small, and the results 
should be interpreted as indications only. 
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TABLE IV 



Ratio of Distribution Graded 
1 for Different Cn Ranges 



Spatial 



Cn Frequency Normal Log Normal 



H 




0 


5/6 


5/6 


I 


1 - 


6 


2/6 


3/6 


G 


7 - 


12 


3/6 


2/6 


H 


13 - 


19 


0/6 


0/6 


M 




0 


3/7 


4/7 


E 










D 


1 - 


6 


1/7 


5/7 


I 










U 


7 - 


12 


0/7 


3/7 


M 


13 - 


19 


0/7 


2/7 






0 


4/7 


4/7 


L 


1 - 


6 


3/7 


5/7 


0 


7 - 


12 


0/7 


1/7 


W 


13 - 


19 


1/7 


2/7 



D. RECOMMENDATIONS FOR FURTHER WORK 

To reduce analysis computing time, the HP -9825 should be 
equipped with additional memory capacity. To eliminate use 
of tape for data storage, memory should be increased to 25K. 

The data analysis for each run should be stored and 
analyzed for mean and standard deviation. Several runs on 
the same data set should be performed to verify distribution. 
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Number of scans sampled in each run should be increased 
by either solving signal trigger problem at low tape speed 
or by reprogramming NIC-80 computer to perform and store 
partial analysis results. 



38 



DATE 31176 TIME 2320 



H 

cn 

PJ 

H 



O' 

on 



SC 

CJ 



oooooooooooooooooooo 

WPJPJPJWPJPJPJPJPJWPJPJPJPJPJPJPJPJtJ 
OH<tr^Lnoocn<j-^Homna\cococNr^o 
^■HOHini^o^vo^-c^Nrr^cn^oinin^-oo 
HOOOCvJHH^O^r^HCOtNOO^Ol^CONfNj- 
MfCMCM^-ICMvOO'^-lMl'UO^-ir'^COvO-Cj _ r^COCOCO 
NOHco^cn^^Lnoo^cMoocnoo^nmncn 

i—l i—l i — I i — I i — I i — I i — I i — I CM i— | CM CM i— | CO i — I i — I i—l i — I CM CO 







i—l 


i—i 


1—1 


1—1 


1 — 1 


I— 1 


1 — 1 


1—1 


1 — 1 


1 — 1 


1 — 1 


1 — 1 


1 — 1 


1 — 1 


1 — 1 


1 — 1 


1 — 1 


1 — 1 


1 — 1 


1—1 






O 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 




<y 

rr\ 


w 


w 


w 


w 


w 


w 


w 


w 


w 


w 


w 


w 


w 


w 


w 


PJ 


w 


w 


w 


w 


CM 


UJ 


CO 


o 


CO 


CO 


o 


CO 


1 — 1 


mT 


VO 


1 — 1 


CO 


co 


mT 


CM 


o 


o 


CO 


00 


o 


CO 




i i 


CM 


co 


CM 


1 — 1 




VO 




uo 


CO 


CO 


co 


uo 


VO 


mT 


<r 


CM 


LO 


o 


CO 


LO 


UO 


r*“i 

nr* 


r^ 


vO 


CM 


LO 


uo 


<f 


vO 


CO 


CO 


CO 


1 — 1 


o 


CM 


O' 


00 


VO 


I—l 


Ml- 


O' 


o 




Pm 

r ■) 






r^ 




co 


CO 






CM 


CM 


r^ 


CM 


CO 


Mf 


00 


00 


uO 




Mj- 


VO 


t — \ 




CO 


1 — 1 


vO 


VO 


VO 


VO 


VO 


CM 


o 


av 


CO 


UO 


CM 


00 


co 


CO 


00 


CO 


r-* 


CM 






i—i 


I—l 


I—l 


1 — 1 


1 — 1 


1 — 1 


1 — 1 


CM 


1 — 1 


1 — 1 


1 — 1 


1 — 1 


1 — 1 


1 — 1 


1—1 


CO 


CM 


1 — 1 


I—l 


CM 



cn 

CJ 

M 

H 

a. 



C 

CJ 



on 

M 

H 



Z 

< 



o 






CM 


CM 


CM 


CM 


CM 


CM 


CM 


CM 


CM 


CM 


CM 


CM 


CM 


CM 


CM 


CM 


CM 


CM 


CM 


CO 




o 




o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 




z 


> 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 






w 


w 


w 


w 


w 


w 


w 


w 


PJ 


PJ 


PJ 


PJ 


PJ 


w 


PJ 


PJ 


pul 


PJ 


w 


w 


PJ 


1 


• 


Q 


CO 


o 


Mf 


1 — 1 


CM 


CO 


LO 


co 


VO 


CO 


r- 


00 


VO 


00 




CM 




CM 


CO 




w 


H 




o 


00 


o 


uo 


o 


o 


uo 


o 


VO 


o 


ON 






CO 


o 


Mj- 


1 — 1 


1 — 1 


uo 


uo 




on 




o 


VO 


VO 


VO 


VO 


Mf 


CM 


CO 


o 


1—1 


o 


CM 


os 


CO 


CO 


1—1 


CO 




uo 


o 




M 


H 


1—1 


r^ 


CO 


CO 


CO 


o 


CO 


VO 


VO 


CM 


o 


1—1 


<r 




1—1 


r- 




1—1 


CM 




040 


Q 


on 


CO 


VO 


CO 


CO 


o 


1 — 1 


ON 


co 


00 


00 


00 


r-» 


VO 




Mj- 


CM 


CM 


1 — 1 


o 


00 


H 




CO 


CM 


CM 


CM 


CM 


CM 


1 — 1 


1 — 1 


1 — 1 


1—1 


1—1 


1—1 


1—1 


1—1 


1—1 


1 — 1 


1 — 1 


1 — 1 


1 — 1 


CO 



1 — 1 


1 — 1 


1—1 


1—1 


1 — 1 


1—1 


1—1 


CM 


CM 


CM 


CM 


CM 


CM 


CM 


CM 


CM 


CM 


CM 


CM 


CM 


o 


o 

1 


o 


o 

1 


o 


o 


o 


o 


o 

1 


o 


o 


o 


o 


c 


o 


o 


o 


o 


O' 


o 


w 


PJ 


PJ 


PJ 


1 

PJ 


PJ 


PJ 


PJ 


1 

PJ 


PJ 


pj 


PJ 


PJ 


pj 


PJ 


PJ 


PJ 


PJ 


PJ 


PJ 


00 


uo 


VO 


Mi- 


Mf 


o 




CM 




LO 


CM 


uo 






CO 


VO 


CO 


1 — 1 


CO 


o 


Mf 


CO 


00 


r-* 


vo 


uo 


o 


CM 


<r 


CM 


uo 


1 — 1 


CO 


o 


VO 


o 


CO 


CO 


VO 


<r 


1—1 


<r 




CO 


CO 


VO 


VO 


o 


CO 


vt 


uo 


CO 


o 


CO 


o 


r- 


CO 


o 


CM 


co 


CM 


CO 


CO 


CO 


CO 


o 


1—1 


sr 


<r 


o 


CM 


1 — 1 


CO 


CM 


Mf- 


CO 


1 — 1 


o 


1 — 1 


CO 






VO 


LO 


CO 


CM 


o 


CO 


uo 


uO 






CO 


Mf 


o 


U0 


<r 


1— 1 


o 


00 



i—l i—l i—l i — I i— I i — I i — I CO P*"* VO LT) Mf CO CO CO CNl CM CM CNl i — I 



OCJ 

H 

< 

PL. 

on 



OHCMroM’invDMoo^OHcvjrONtmvorNOO!^ 



39 



25991 



MERN HNP 4- — I STHNPRRI> 



Lrl 

Z - 

□ - 



h 

H 



Ifl 

z 

- E 

> V 

u in 

Cx 

L 

□ 

□ 

Z 



□ 

N 

PI 

N 

hi 

X 



L 

h 

L 



in 

t" 



FI 

U 

I- 

E 

Q 



JM 

h 

U 

r 

□ 

in 

v 

h 

L 

□ 



□ 

N 
PI 

in 

Q 

X 
1- 

in a 

Z T 



U 

J 

u 

> 

E 

Z 



l 

U 



a 



z 

K> 



\A 



a 



\A 



a 



amnihini/irFiN-H 

— SSQSSQSSSS 



JIN 



» » 



40 



SPRTIRL FREQUENCY IN 1/ MILLIRRPIHNE 
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UNITE - 1/ 3 ETRNI>RRD OEVIflTIDNS 



APPENDIX A 



Operating Instructions for Data Analysis 
A. Initialization 

1. NIC-80, Analysis program not loaded 

(a) Turn on computer with key, depress "STOP" 

(b) Turn on Data Bus unit 

(c) Set toggle switches 0 - 11 to up position 

(d) Depress "LOAD PC" 

(e) Depress "CONTINUE" 

(f) Press "EXECUTE" 

2. HP-9825 

(a) Turn on HP -9 825 

(b) Insert Data Analysis cassette and load track 0, 
file 1 for NIC-80 loading program 

(c) Push "RUN"; c, d, i, r, n, o, 1, s, t is displayed. 

(d) Push "n" and "CONTINUE"; load file # is displayed. 
The NIC-80 data analysis program is located on files 2, 3, 4 
and 5. Enter one of these numbers and push "CONTINUE". 

(e) Display - c, d, i, r, n, o, 1, s, t. Push "1" 
and continue. The section of the program will be transferred 
to NIC-80. At the end of each transfer, check sums are 
printed and should be identical. 

(f) Repeat steps (d) and (e) until entire program is 
transferred. If NIC-80 contained MTF analysis program pre- 
viously, only files 2 and 3 are necessary to transfer. 
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(g) On NIC-80, push "STOP", depress "START" and push 
"EXECUTE" on completion of program loading. 

3 . Tape Recorder 

(a) Turn switch to 3.75 IPS 

(b) Switch on DC bias boxes. These should be turned 
off during non-operating periods to avoid draining batteries . 

(c) Set frequency control knob on appropriate channels 
to FM and lkc bandpass for signal and trigger. Voice channel 
requires 37.5 IPS and 100 Kc for operation. 

(d) Turn on oscilloscope. 

4. Pulse Generator 

(a) Turn on power 

(b) Set to external trigger pulse. 

(c) Set width to 10 microseconds. On tape recorder, 
depress "DRIVE" and adjust oscilloscope for optimum display. 

B . EXECUTION 
1. HP-9825 

(a) Load track 0, file 8 and push "RUN" 

(b) Data entry; 

Wavelength, e.g., .6328 
Optics , e.g. , .4572 
Dwelltime, e.g., 95 
Number of scans, e.g., 125 

Push "CONTINUE" after each entry. On last "CONTINUE" NIC-80 
will enter Display/ Standby mode. HP-9825 will display 
"CONTINUE WHEN READY". 
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2. NIC-80 

(a) Verify green control light is on 

(b) Adjust CRT for optimum display of scans 

(c) Adjust DC level by knobs on bias boxes. Use low 
resolution on vertical display scale, and adjust signal to 
center around horizontal scale. 

(d) Adjust horizontal position of signal by pulse 
generator delay knob so as to center around vertical axis . 

3. To start data collection, push "CONTINUE" on HP-9825. 
When all scans are completed, NIC-80 will return to Display/ 
Standby mode and may be taken off-line. 

4. Data Analysis (HP-9825 only) 

(a) On completion of data collection, "HISTOGRAM AND 
CONTINUE" is displayed. Check that plotter and printer are 
ready before pushing "CONTINUE". 

(b) Data entry; 

Date, e.g., 021577 (month, day, year) 

Time, e.g. , 1423 
Cn arg, e.g. , 8.07 
Power, e.g., -8 

Push "CONTINUE" after each entry. 

(c) Unnormalized values for MTF spatial frequencies 0 
and 1 are printed for verification of absolute and relative 
values . 

(d) Data entry; "No. of class intervals" is displayed. 
Enter number of histogram bins per unit standard deviation. 
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Normally 3. Push "CONTINUE”. 

(e) Display "FREQ PLOT, YES-1". Push "1" and "CON- 
TINUE" if plot of MTF mean + 1 standard deivation for spatial 
frequency 0 - 19 is required. 

(f) Display "LN DATA? YES-1". Enter 1 and push 
"CONTINUE" if data analysis of normal distribution is completed. 

Display "LN ANALYSIS". Push "CONTINUE" if data 
is to be converted to Ln data. 

(g) Display "DATA SUM? YES-1". Enter 1 and push 
"CONTINUE" if data summary is required. 

Display "SUMMARY, CONT.". Push "CONTINUE" to 
start execution of data summary program. 

(h) Display "SPAT FREQ". Enter spatial frequency 
number (0 - 19) for which histogram will be plotted. 

Display "REPLACE PAPER AND CONTINUE". Push 
"CONTINUE" when plotter is ready. 
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0: if X=l;gto "START" 

1: mdecjpar l;frat l,fzl0.5;enp "wavel " , W, "opt " ,0 
2: enp " DWELLTIME" , rO , " # OF SCANS",N 

3: 20-P;dim F [P] , A [P] ,M [P] ,S [P] ,X [P] ; wr t 2 .1, W,0, rO ,N; oar 0 
4: for J=1 to 20 ; rdb ( 2) -A; rdb (2 ) -3 ; 512* B+A-F [ J] 

5: if F[J]>524287 ;F[J]-108576-F[J] 

6 : next J ; tr k 1 

7: enp "CONTINUE WHEN READY", rO 
8 : wtb 2 , r 0 

9: for 1=1 to N;for J=1 to 20 ;rdb (2 ) -A; rdb (2 ) -3; 512*B+A-A [ J] 
10: if A[J]>524237 ;A[J] -104S575-A[ J] 

11: rdb (2) -A; rdb ( 2) -B; 512*3+A-r0 
12: if r0>524287 ; r0-1048576-r0 
13: A[J] *A[J] +r0*r0-r0 ;/r0-A [J] 

14: A[J]/F[J]+A[J] ;A[J]*A[J]+S[J]+S[J] ;M [ J ] +A [ J ] -M [ J] 

15: next J;rcf I,A[*] 

16 : next I 

17: dsp "HISTOGRAM ?,and cont i nue " ; stp ;flt 5 

13: enp " DATE" ,r 99 , "TIME " , r 98 , "CN ARG" , r9 4 , " POWER" , U 

19: "PLOT": for Q=1 to 20 ; tr k 1 

20: (S[Q]-M[Q]*M[Q]/N)/(N-1)-S[Q] ?M[Q] /N*M [Q ] ;/S[Q]*5[Q] 

21: next Q 

22: fit 5 ; pr t "MTF 0" , M [ 1 ] , "MTFl " , M [ 2 ] ; if r97 = l;gto "DS" 

23: enp "NO CLASS INTERV" , T ; beep 

24: "START" : enp "FREQ PLOT , YES = ’l" , r9 2 ; if r92 = l;gto " FREPLO" 
25: if r97=l;gto "DS" 

26: "L" : enp " LN DATA? , YES = 1 " , r 97 ; i f r 9 7 = 1 ; 1-Z ; g to "LN" 

27: "DS" :0-X;enp "DATA SUM? , YES= 1" , X ; if X=l;tr!< 0;ldf 9 

28: "PLOTO" : enp" "SPAT FREQ" , A ; A+ 1-Q ; tr k 1 

29: dsp "REPLACE PLOT PAPER AND CONTINUE beep ; stp 

30 : M [Q] -3*S [Q]-C;M [Q] +3*S [Q] *D 

31: S [Q] /T-R;C-2 ,2*R-r91 ;cs iz 1.6, 1.5,1 

32: scl C-2*R, D, - .2 , 1 .7 ; pi t D,0,l;plt C-R/10 , 0 , 2 ; f xd 1 

33: pit r91+R/2,0,l; pit r9 l + R/2 , . 7 , 2 ; f or 1=0 to .7 by .1 

34: ltr r91+R/l. 5, I; lbl I;plt r9 1 + R/2-R/10 , 1 , 1 

35: pit r91+R/2+R/10 , 1 , 2 ; next I 

36: -3*T-r 88; fxd 0 

37: for I=C to D by R;olt I,.01,l;plt I, -.01, 2 

38: ltr I-R/3 , - . 03 ; lbl" r9 8 

39: r88+l-^r88; next I;pen;csiz 1.6, 1.5,1 

40: ltr C+3*R, 1.6 ; lbl "RELATIVE DISTRIBUTION OF MTF" 



Figure 8A. Data Analysis Program Listing 
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41: pen; pit M[Q],0;plt M[Q],,7,2 

42: ltr C+R,-.15;lbl "UNITS = 1/",T," STANDARD DEVIATIONS" 

43: "CURVE" :for I=C to D by ( D-C) /100 ; ( I-M [Q] ) /S [Q]+r 1 

44: rl*rl/2+rl; exp (-r 1) *R/ (S [Q]*{ ( 2 * tt ) )*r2 

45: pit I,r2,-2;next I ;C*r 2; r2-R-*-r 3 

46: C+r 2 ; r 2 + R+r 3 ; f It 5;A-*-A 

47: for K=ll to 1 1+T*6 ; 0-*-r K; nex t K 

43: "HIST" : for 1=1 to N;C-r2 ; r2+R-r3 

49: ldf I , A [ * ] ; for K=ll to 11+T*6 

50: fit 5 ; if A[Q]>r2;if A [ Q] < r3 ;r K+l-r K 

51: r3-*-r 2 ; r3+R-»-r 3 ; nex t K;next I 

52: C+r 2 ; r 2+R+r 3 ; f or K=ll to 11+T*6 ; r K-*r 1 

53: r l/N-*-r l;pen;plt r2,0,-2;plt r2,rl,2;plt r3,rl,2 

54 : pit r3, 0, 2 

55 : ( r 2-M [ Q ] ) / S [Q]-r7; rl*rl/2*rl ; exp ( -r 7 ) *R/ ( S [ Q] * V ( 2 * tt ) ) -r4 

56 : (r 3-M[Q] )/S[Q]*r7; r7*r7/2-r7 ;exp (-r7) *R/(S [Q]*V ( 2 * tt ) ) + r5 

57 : R*r 4 +abs (r5-r4)/2*R-*-r3;r8/R+r8 

58 : (N*r8-rl*N)~2/(N*r8)-*-r9;r9 + rlO-*-rlO 

59: r 3-*-r 2 ; r3 + R-*-r 3 ; next K 

60: "LABEL" :C + R/3-rl; fxd 0;ltr rl,l,45 
61: lbl "SPATIAL FREQUENCY" , Q-l 

62: if r97=l; 1 tr rl, 1.25; lbl "TEST DIST LN NORMAL" ;gto +2 
53: ltr rl, 1.25; lbl "TEST DIST NORMAL" 

64: ltr rl, 1.35; lbl "NO OF SCANS", N 

55: fxd 4; ltr rl, 1.15; lbl "NAVE LENGTH " , W; 1 tr rl,1.05 
66: lbl "OPTICS" ,0 

67: fxd 0 ; ltr rl,.9;lbl " DATE " , r99 , " TIME" , r 98 ; f It 5 
53: M[Q] +R/3+rl; ltr rl,1.35;lbl "St. dev",S[Q] 

69: ltr rl, 1.45; lbl "MEAN" , M [ Q] 

70: ltr r 1, 1 . 25 ; lbl "CHI SQ",rl0;ltr rl,1.15;fxd 3 
71 : lbl "CN" ,r 94 ; f xd 0 
72: lbl " E " , U 

73: for 1=1 to 90 ; 0+ r I ; nex t I 
74: if r97=l ; g to "DS" 

75 : gto "L" 

76: "LN" :dsp "LN At4ALYSIS " ; stp ;trk 1 

77: for J=1 to 20; 0+M ( J] +3 [ J] ;next J 
7 3 : for 1=1 to N 

79: ldf I,A[*];for J=1 to 20;if A[J]=0;if J=20;gto +5 
80: if A[J]=0;next J 



Figure 8B. Data Analysis Program Listing 
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31: ln(A[J] )*A[J] 

32: A [J]*A[J] +S [J] -S [J] ; M ( J ] +A [ J] -M [ J ] 

33: next J;rcf I,A[*] 

34: next Ijgto "PLOT" 

35: "FREPLO" : scl -3 , 20 , - . 2 5 , 1 . 5 ? f xd 0;csiz 1.6, 1.5,1 
86: pit 20 , 0 , 1 ; pi t 0,0,2;plt -1.35 ,0 ,1 ;plt -1.85,1,2 
37: for 1=1 to 20; pit 1-1 , . 01 , 1 ;o It 1-1, - . 01 , 2 ; nex t I 
38: for 1=1 to 20 by 5 ; I-l-r95 ; ltr r 95- . 3 , - .0 5; lbl r95 
8 9 : next I ; f xd 1 

90: for 1=0 to 1 by ,1,-ltr -1.35,I;lbl I;plt -1.9, 1,1 

91: for 1=0 to 1 by ,l;ltr -1.35,I;lbl Ijplt -1.9, 1,1 

92: pit -1 ,8 , 1 , 2 ; next I 

93: ltr 2, -.19; lbl "SPATIAL FREQUENCY IN 1/ MILL I RAD IANS " 
94: ltr 2, 1.45; lbl "MTF MEAN AND + - 1 STANDARD DEVIATION" 
95: fxd 0;ltr 2,1.35 

96: lbl " DATE " , r 99 , " TIME",r98," MO OF SCANS", N 

97: fxd 4; ltr 5, 1.25; lbl "WAVELENGTH", W, " OPTICS", 0 
93: fit 5 ; 1 tr 5,l,15;fxd 3;lbl "CN",r94;fxd 0;lbl "E",U 
99: M [ 1 ] ♦ r 9 3 

100: csiz 3 , 1 , 1, 90 ; ltr -2.1,.4;lbl "MTF" 

101: for Q=1 to 20 ; 50 +Q*r 43 ; 7 O + Q* r 49 ; M [Q ] /r 93-*- r 43 
102: S [Q ] /r 93*r 49 

103: pit Q-l, r 43-r49 , 1; pit Q-l , r48+r49 , 2 

104: pit Q-l ,05 , r 43 , 1 ; pi t Q- . 9 5 , r4 3 , 2 ; pen ; nex t Q;gto "L" 



Figure 8C . Data Analysis Program Listing 
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0: asp "SUMMARY ,CONT" ; beep ;stp 
1: "X" : f or Q=1 to 20 
2 : S [Q] /T-R 

3: A [Q]-3*S [Q]-C; M[Q] +3*S [Q] -D 
4: for K= 11 to 11+T* 6 ; 0 + r K ; next K 
5: "HIST":for 1 = 1 to N; C*r2 ; r2+R-*r 3 
6 : tr k 1 

7: ldf I , A [ * ] ; f or K=ll to 11+T*6 

8: fit 5; if A[Q]>r2;if A [ Q ] < r 3 ; r K+ l-»r K 

9: r 3->r 2 ; r 3 + R*r 3 ; nex t K;next I 

10: Or2; r2+R*r 3 ; f or K=ll to 11+T*6 ; r K-*r 1 

11: rl/N-*-rl 

12 : (r2-M [ Q ] ) /S [Q]*r 7 ; r7 *r7/2-r7 ; exo (-r7 ) *R/(3 [Q] * / ( 2*ir ) ) *r4 
13: (r3-M[Q])/S[Q]*r7;r7*r7/2*r7;exp(-r7)*R/(S[Q]*V(2*ir))-r5 
14 : R*r4+abs(r5-r4)/2*R-r3; r3/R-r3 
15: (N*r3-rl*N)~2/(N*r3)-*-r9;r9 + rl0-»-rl0 

16: r3+r 2; r3 + R-*-r 3 ; next K ; rlO*X [Q] ;0+r 10 ; next Q;gto "PRINT" 

17 : "PRINT" : fmt 2, 25x , "DATE" , 2x , f 6 ,0,2x, "TIME" , 2x, f4 . 0 

13 : fmt 1, 2 5x , "Cn" , 2 x , f 5 , 3 , 2x , "E " , f 3 . 0 , 2x , "OPTICS" , 2x , f 6 .4 , 2x 

19 : fmt 3 , 3 5 x , " DATA SUMMARY " 

20: if 2=1 ; g to "LI" 

21: wr t 6 .3, 13 
22: wtb 6,10,10 
23: wr t 6,2,r99,r98 
24: wtb 6,10,10 
25: wrt 5 . 1 , r 94 , 0 

25 : wtb 5,10,10 

27: fmt 4 , 2 5x , " TEST D 1ST .NORMAL" , 2 6 X , " TEST DIST LN NORMAL" 

28 : wr t 6.4 
29 : wtb 6 , 10 , 1 G 

30: fmt 5 , 1 5 x , " S P AT FR" , 5 x , "MEAN " , 1 lx , " ST . D EV" , lOx , " CH I 3Q",10x,"CHI 
31: wrt 6,5;wtb 6,10,10 

32: fmt 5 , 15x , f 2 . 0 , 6 x , el 1 . 5 , 5 x , e 1 1 , 5 , 5x , e 1 1 . 5 

33: for 1=1 to 2Q;wrt 6 .5 , I- 1 , M [ I ] , S [ I ] , X [ I ] ; nex t I 

34: for 1=1 to 20; wtb 6, 27, 10; next I 

35: fmt 7,72x,ell.5 

36 : if 2=1; gto "LI" 

37 : tr k 0 ; ldf 8 

38 : "LI" : fmt 7 ,72x,ell .5 

39: for 1 = 1 to 20;wrt 6 . 7 ,X [ I ] ; nex t I 

40: tr k 0 ; ldf 8 



Figure 9. Data Summary Program Listing 
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0 


110600 


MEMA (600 /READOUT LIGHTS 


1 


4304 


LCWORD /CONTROL WORD 


2 


44034 


STU /READ STATUS 


3 


10377 


ANDA (377 /COMPARE SETTINGS 


4 


462340 


A— MZ (340 /CHECK HARD SETTING 


5 


0 


JMP 0 /RETURN CHECK LOOP 


6 


2175677 


ZERMA 1677 /0 TO ACC & 1677 


7 


4304 


LCW (0 /LOAD CONTROL WORD 


10 


3001737 


JMS@ INP /INPUT 


11 


3001745 


JMS@ 1745 /PUT AC 


12 


1616 


1616 /ADR OF WAVELENGTH 


13 


3001737 


JMS@ INP /INPUT 


14 


3001745 


JMS9 1745 /PUTAC 


15 


1614 


1614 /ADR OF OPTICS 


16 


3001737 


JMSI§INP /INPUT 


17 


3001745 


JMS91745 /PUTAC 


20 


1600 


1600 /ADR OF DWELLTIME 


21 


3001737 


JMSdINP /INPUT 


22 


3001745 


JMS81745 /PUTAC 


23 


1620 


1620 /ADR NO OF SCANS 


24 


110020 


MEMA (20 /COUNT VALUE 


25 


2405677 


AC CM 1677 /COUNT INDEX 


26 


2001500 


JMS 1500 /CAL ADRESS 


27 


3111576 


MEMA@ 1576 /CAL VALUE 


30 


10777 


ANDA (777 /9 BITS 


31 


6474 


HSP /IS HP 9825 READY 


32 


31 


JMP 31 /WAIT LOOP 


33 


4473 


PHS /WRITE TO HP 9825 


34 


3111576 


MEMA@ 1576 /CAL VAL RECALLED 


35 


405031 


RIS 11 /RIGHT SHIFT 


36 


64 74 


HSP /IS HP 9825 READY 


37 


36 


JMP 36 /WAIT LOOP 


40 


4473 


PHS /WRITE REMAINING DATA 


41 


2125576 


MPOM 1575 /INCR. CAL ADRESS 


42 


2703677 


MMOZ 1677 / DECR . & TEST 


43 


27 


JMP 27 /RETURN LOOP 


44 


3001742 


JMS@1742 /GETAC 


45 


1600 


1600 /ADR. OF DWELLTIME 


46 


3001715 


JMS? 1715 /FIX 



Figure 10A. Nic 80 Program Listing 
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47 

50 

51 

52 

53 

54 

5 5 

56 

57 

60 

61 

62 

63 

64 

6 5 

66 

67 

70 

71 

72 

73 

74 

75 

76 

77 

100 

101 

102 

103 

104 

105 

106 

107 

110 

111 

112 

113 

114 



3111722 

4311 
110000 

4312 
2111711 

4304 

4302 

5220 

6464 

53 

44463 

462314 

64 

3001723 
2111704 
2405677 
2111626 
2405672 
2111712 
4304 
4302 
52 20 
100 
0 
0 

110200 

4304 

2111700 

2405677 

2711701 

2405676 

110040 

2405675 

2165674 

3111677 

2505674 

3111676 

2505674 



MEMA@1722 / FACM 

LDWELL /LOAD DWELL TIME 

MEMA ( 0 /0 TO ACC 

LDELAY /LOAD DELAY TIME 

MEMA 1711 /STAND BY LOOP 

LCWORD /LOAD CONTROL WORD 

SET /SET MEASURE FF 

STOP /READ 1 SCAN 

HSRF /SKIP ON HP REQUEST 

JMP 53 /RETURN STAND BY LOOP 

RHSR /READ HP 

A-MZ ( 314 /L FOR FUTURE BRANCH 
JMP 64 /FUTURE BR. ADR. 

JMS@ 172 3 /SHORT TERM RESET 
MEMA 1704 /B4 ( 0 ) 

ACCM 1677 / 

MEMA 1626 /# OF SCANS 

ACCM 1672 / 

MEMA 1712 /CONTROL WORD #2 
LCWORD /LOAD CONT.WORD 
SETM /SET MEASURE FF 
STOP /READ 1 SCAN 
JMP 100 / 

0 /NOT USED 
0 /NOT USED 
MEMA (200 / 

LCWORD /DISPLAY CONTROL 
MEMA 1700 /30 (0), LEFT END 

ACCM 1677 / 

MMOA 1701 /B0 (1777) .RIGHT END 

ACCM 1676 / 

MEMA (40 /32 POINTS 



ACCM 1675 
ZERM 1674 
MEMA @16 77 
A+MM 1674 
MEMA01676 
A+MM 1674 



/ 

/ 

/INTEGRAL 

/ 

/ 

/ 



BASE LINE 



Figure 10B. Nic 80 Program Listing 
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115 

116 

117 

120 

121 

122 

123 

124 

125 

126 

127 

130 

131 

13 2 

13 3 

134 

135 

136 

137 

140 

141 

142 

143 

144 

145 

146 

147 

150 

151 

152 

153 

154 

155 



2125677 


MPOM 1677 / 


2705676 


MMOM 1676 / 


2707675 


MMOMZ 1675 / 


111 


JMP 111 / 


2111674 


MEMA 1674 / 


5026 


RASH6 /DIVIDE BY 6 


2405674 


ACCM 1674 /AVERAGE DC COMPONENT 


2111700 


MEMA 1700 /B0 (0 ) 


2405677 


ACCM 1677 / 


3111677 


MEMA01677 / 


2471674 


A-MA 1674 /SUBTRACTING 


3405677 


ACCM@16 77 / 


2135677 


MPOM A 1677 / 


2463701 


A-MZ 1701 / 


126 


JMP 126 /RETURN ADR. LOOP 


110002 


MEMA ( 2 


3001722 


JMSQ1722 /MOVE B0*B2 


2111703 


MEMA 1703 /B3 ZERO LOOP 


2405677 


ACCM 1677 / 


2111700 


MEMA 1700 /B0 


2405676 


ACCM 1676 / 


3165676 


Z ERM@ 16 76 / 


3165677 


ZERM@1677 / 


2125677 


MPOM 1677 / 


2135676 


MPOM A 1676 / 


2463701 


A-MZ 1701 / 


142 


JMP 142 /RETURN ADR. LOOP 


3001776 


JMS01776 / 


3001622 


JMS@16 22 / DFT DATA TRANSFER 


4032 


ASCNTR /ADVANCE SWEEP COUNTER 


2703626 


MMOZ 1625 /DECREASE SCAN COUNT 


71 


JMP 71 /RETURN ADR. NEXT SCAN 


53 


JMP 53 /RETURN ADR. STAND BY L. 



Figure 10C. Nic 80 Program Listing 
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3445 


15 2 


152 /RETURN ACRES S 


3446 


110024 


MEMA ( 24 


/# OF SPAT . F R . 


3447 


2405675 


AC CM 1675 


/COlNTING INDEX 


3450 


2111702 


MEMA 1702 


/B2, REALS 


3451 


2405677 


ACCM 1677 


/CURRENT B2 


3452 


2111703 


MEMA 1703 


/B3, IMAGINARIES 


3453 


2405676 


ACCM 1676 


/CURRENT B 3 


3454 


3111677 


MEMA01677 


/VALUE B 2 


3455 


10777 


ANDA (7 7 7 


/9 3 ITS 


3456 


64 74 


HSP /HP 


READY 


3457 


1456 


JMP 1456 


/WAIT LOOP 


3460 


4473 


PHS /WRITE TO HP 


3461 


3111677 


MEMA01677 


/RECALL B2 


3462 


405031 


RIS 11 /RIGHT SHIFT 


3463 


6474 


HSP / 




3464 


1463 


JMP 1463 


/ 


3465 


4473 


PHS / 




3466 


3111676 


MEMA91676 


/ 


3467 


10777 


ANDA (777 


/ 


3470 


6474 


HSP / 




3471 


1470 


JMP 1470 


/ 


3472 


4473 


PHS / 




3473 


3111676 


MEM A@ 16 7 6 


/ 


3474 


405031 


RIS 11 / 




3475 


6474 


HSP / 




3476 


1475 


JMP 1475 


/ 


3477 


4473 


PHS / 




3500 


2125677 


MPOM 1677 


/INCREASE 32 


3501 


2125676 


MPOM 1676 


/INCREASE 33 


3502 


2707675 


MMOMZ 167 


5 /DECREASE & TE 


3503 


14 54 


14 54 /RETURN ADR. 


3504 


1001445 


JMP @144 5 


/RETURN TO MAIN 



Figure 10D. 
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